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ABSTRACT 


The  heeling  behavior  of  a  fully-submerged  foil 
hydrofoil  craft  turning,  has  been  studied  by  producing  a 
digital  ccmputer  simulation  program,  using  a  Digital 
Simulation  Language  on  the  digital  computer  IBM  system 
360/67. The  High  Point  PC(H)-1  is  used  as  the  model  of  this 
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!•   INTRODUCTION 


This  thesis  deals  with  the  yaw  motion  and  turning 
characteristics,  in  calm  water,  of  the  hydrofoil  craft 
PC(H)-1  High  Point.  This  will  be  the  model  used  in  the 
computer  simulation  program  for  this  study. 

Iu  previous  theses,  the  system  was  made  conditionally 
stable  in  two  modes,  the  pitch-heave  mode  and  the  roll  mode. 
In  this  thesis,  keeping  the  system  stable  and  working  in  six 
degrees  of  freedom,  its  behavior  will  be  analyzed  in  three 
modes,  pitch  mode,  roll  mode  and  yaw  mode,  when  the 
simulated  model  is  turning  by  the  effects  of  a  helm  step 
command,  acting  on  the  astern  flaps.  This  analysis  will 
also  include  the  study  of  rudder  effects,  aft  starboard  and 
aft  port  flaps  effects  and  the  heeling  behavior  of  the 
system.  In  a  further  investigation  the  forward  flap  will  be 
split  in  two  parts  and  pitching  and  rolling  effects  will  be 
analyzed. 

Since  slight  changes  in  the  manual  control  of  the 
system  occur  in  such  short  time  intervals  and  provoke 
tremendous  effects,  the  system  necessitates  an  automatic 
control  that  will  provide  the  necessary  outputs  to  operate 
the  model  according  to  the  given  inputs.  This  automatic 
control,  already  existent,  will  be  improved  and  adapted  to 
the  new  motion  of  the  craft. 

Finally  the  favourable  roll  angle  for  a  certain  speed 
and  a  certain  helm  step  command  will  be  calculated  and 
compared  with  the  actual  roll  angle  of  the  craft. 


II*   HYDROFOIL  EQUATIONS  STUDY 


A.   FOIL  TYPES  AND  FOIL  ARRANGEMENT 

There  are  two  principal  foil  types  used  in  the 
construction  of  hydrofoil  crafts,  the  surface-piercing  foil 
and  the  fully-submerged  foil, as  shown  in  Figure  1. 


Surface-Piercing  Foil  Fully-Submerged  Foil 

FIGURE  1 
TYPE  OF  HYDROFOIL  CRAFT 
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The  surface-piercing  foil  provides  a  self-stabilization 
for  the  craft  in  pitch,  height  and  roll.  On  the  other  hand 
the  fully-submerged  foil  cannot  adjust  itself  to  a  changing 
lift.  This  can  only  be  achieved  by  a  change  in  incidence 
angle  of  the  foil  flaps.  The  fully-submerged  foil  requires  a 
more  sophisticated  sensing  arrangement  to  control  that 
incidence. But  the  fully-submerged  foil  has  a  great 
advantage,  when  the  craft  is  operating  in  heavy  sea. The 
surface-piercing  foil  responds  to  relatively  small  changes 
in  the  surface  of  the  sea  and  gives  a  very  rough  ride. On  the 
other  hand,  the  fully-submerged  foil  is  much  less  affected 
by  the  effects  of  the  surface  of  the  sea. This  paper  will 
only  deal  with  the  fully-submerged  foil  craft. 

The  model  for  this  simulation,  the  High  Point 
PC(H)-1,uses  the  foils  in  canard  arrangement,  which  has  a 
secondary  foil  as  forward  foil  and  a  main  foil  as  aft  foil. 
Both  foils  are  supported  by  vertical  struts.  Figure  2, 
important  for  the  understanding  of  the  digital  computer 
simulation  program,  shows  the  configuration  of  the  forward 
and  aft  foil. The  drawing  and  the  dimensions  of  the  craft  and 
the  foils  are  shown  in  Figures  3A  and  3B. 
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FOILS  CONFIGURATION 
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B. 


COORDINATE  SYSTEMS 


Three  coordinate  systems  will  be  used  in  this  study  of 
the  hydrofoil  craft  dynamics.  They  are  :  the  Body  Axes,  the 
Earth  Axes  and  the  Water  Axes.  Each  of  these  systems  is  a 
right  hand  orthogonal  system. 

The  body  axes  coordinate  system  (X,Y,Z),  shown  in 
Figure  4,  has  its  origin  at  the  center  of  gravity  of  the 
craft  and  is  fixed  relative  to  the  craft  with  the  X  axis 
running  longitudinally  through  the  craft.  This  system  is 
initially  coincident  with  the  origin  of  the  earth  axes 
system,  but  it  moves  with  the  craft  as  time  progresses.  X  is 
positive  forward,  Y  is  positive  to  starboard  and  Z  is 
positive  downward.  ■ 


CrsKwent  BOYf) 


X  -  AXIS 


■*T£D.) 


FIGURE      4 

ORIENTATION  OF  THE  BODY  AXES  WITH  RESPECT  TO  THE 
CRAFT  AND  DIRECTIONS  OF  THE  POSITIVE  VELOCITIES 
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The  earth  axes  coordinate  system  (X  ,  Y  ,  Z  )  is   fixed 

E    E    E 

relative   to   the   earth's  surface.  The  origin  may  be  at  any 

convenient  location,  provided   the   body   system   origin   is 

initially   at   the   same   point.   The   X    axis   lies  in  the 

E 

horizontal  plane  ;  it   is   initially   coincident   with,   and 

positive   in   the  same  direction  as,  the  body  X  axis.  The  Y 

E 

lies  in  the  horizontal  plane  and  is   positive   to   starboard 

when   the  observer  is  facing  the  positive  X   direction.   The 

E 

Z   axis  is  normal  to  the  horizontal  plane   and   is   positive 
E 

downward.   The   craft's   motion   relative   to  the  calm  water 

surface,  (constraint  of  this  paper) ,  is  described   in   terms 

of  this  coordinate  system. 

In  general,  the  body  axes  are  displaced  from  the  fixed 
earth  axes  by  the  three  Euler  angles,  G  (pitch  angle) ,  $ 
(roll  angle)  and  ^  (yaw  angle) .  Figure  6  defines  the 
positive  directions  of  the  Euler  angles.  The  transformations 
between  earth  and  body  axes  are  outlined  in  Ref.3,  by 
transformation  matrices,  however,  only  the  results  will  be 
considered  here. 

The  water  axes  coordinate  system  is  aligned  with  the 
relative  water  velocity  vector  and  resolves  into  lift,  drag 
and  side  force  directions.  All  hydrodynamic  model  test  data 
is  taken  and  presented  as  plots  in  terms  of  this  coordinate 
system. 

The  hydrodynamic  force  and  moment  data  are  obtained  in 
a  water  axis  system  which  is  always  oriented  with  respect  to 
the  relative  water  velocity.  These  data  must  be  transformed 
into  the  body  axes  in  order  that  their  effect  will  be 
properly  included  in  the  equations  of  craft  motions. 

In  general,  the  water  axis  orientation  is  always 
changing  with  respect  to  the  body  axes.  This  leads  to  the 
need  for  a  transformation  matrix,  also  outlined  in  Ref.3,  to 
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resolve  tbe  lift,  drag  and  side  force  values  in  the  body 
axes.  The  orthogonal  water  axis  coordinate  system  is  defined 
in  Figure  5.  The  angles  of  rotation  between  water  and  body 
axes  are  a  and  (3     and  will  be  defined  later. 


vrr 


•o—D 
DfiAd 


FIGURE   5 
ORIENTATION  OF  WATER  AXES  COORDINATE  SYSTEM 
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FIGURE  G 
DEFINITION  OF  POSITIVE  DIRECTIONS  OF  EULER  ANGLE 
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C.   EQUATIONS  OF  MOTION 

The  development  of  the  six  equations  of  motion  for  a 
hydrofoil  follows  the  same  procedure  as  those  for  a 
displacement  type  hull.  This  development  is  clearly  outlined 
in  Ref . 3  and  therefore  will  not  be  done  here. 

The  equations  of  motion  are  summarized  in  this  form  : 

Translation 


U  =  (1/m)  F   -  QW  +  EV  (2-1) 

A 


V  =  (1/m) F   +  PW  -  RU  (2-2) 


H  =  (1/m)  F   +  QU  -  PV  (2-3) 

Z 


Rotation 


P  =  1/1    [l-QR(I   -I    )+(E  +  QP)I   1  (2-4) 

XX  L       ZZ  YY  XZj 

r  22-1 

Q  =  1/1    M-RP(I   -I    )-(P  -R  )I  (2-5) 

YY  L       XX  ZZ  XZJ 

R  =  1/1    |N-PQ(I   -I    )-(QR-P)I 

ZZ  L       YY  XX  XZJ 


(2-6) 


Investigations  with  the  craft  simulation  show  that  the 

terms  RV,PV,  and  PW  are  very  small  and  may  be  neglected. 

Similarly,  the  term  I    is  less  than  10%  of  I    and  less 

XZ  XX 

than  U%    of  I   .   Consequently,  the  terms  containing  I    may 
ZZ  XZ 

be  neglected  for  this  craft. 

Equations  (2-1)  through  (2-6)  are  linear  and  rotational 
accelerations  in  body  axes. 

Given  the  previous  equations,  the  velocity  terms  in 
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body  axes  are  obtained  by  integrating  each  equation  one 
time. 


U  =  J     6 


dt  (2-7) 


V  =  /  V  dt  (2-8) 

H  =  /  H  dt  (2-9) 

P  =  /  P  dt  (2-10) 

Q  =  J    Q    at  (2-11) 

R  =  /  R  dt  (2-12) 

In  order  to  calculate  the  craft's  position  in  earth 
axes  coordinates,  these  velocity  outputs  must  be  transformed 
from  body  axes  coordinates  into  earth  axes  coordinates.  The 
transformation  of  the  linear  equations  become  : 

U      =   Ucos0cos»/AfV  (cosi/'sintfsiiKj!'  -sin  \p  cos<t>) 
E 

+  W  (cos*/'sin0cos<£+sini/'sin  4>  )  (2-13) 


V   =  Ucos^sin^+V (cos^cos^+sin ^  sin0sin#) 
E 

+W  (sinypsindcosQ-cos  yp  sin<t>)  (2-14) 

W      =   -Usin0-«-Vcos0sin4>+Wcos0cos  4>  (2-15) 

E 


Only  the  craft's  velocity  in  the  downward  direction,Z  , 

E 

is  required  to  calculate  its  height  above  the  water  surface. 
That  is,  the  vertical  position  in  earth  axes  follows: 


Z   =  /  W    dt  (2-16) 

E    J         E 
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The  craft's  position  in  earth  axes  coordinates  is 
obtained  from  : 


X   =  /  U   dt  (2-17) 

2    J        E 


Y   =  /  V   dt  (2-18) 

E    J         E 


The  Euler  angles,  4>  ,  0  and  4*  are  also  of  interest 

« 
and  can  be  found  by  integration  of  Euler  rates  <t>  ,  0   and  ^  . 

The  rates  in  turn  must  be  derived  from  the  body  angular 

rates.  The  Euler  rates  are  not  easily  obtained  because  they 

occur  about  axes  which  are  not  orthogonal.  This  fact  can  be 

appreciated  by  referring  to  Figure  6  and  recalling  how  the 

angles  were  defined,  t  is  a  rotation  about  the  Z  axis,  6 

is  a  rotation  about  the  I  axis,  and  4>  is  a  rotation  about 

the  X  axis.  They  result  in  the  following  set  of  equations  : 

<t>  =   P+^sin  6  (2-19) 

6  =    Qcos4>~Rsin  <f>  (2-20) 

■i  -  (Qsin4>+Rcos<2>)  cos0+  (0-P)sin0  (2-21) 

The  Euler  angles  are  obtained  by  integration  of  these 
equations. 

Studies  of  this  craft  indicated  that  during  abnormal 
situations,  roll  angles  of  15  degrees  and  larger  may  occur, 
but  under  normal  operating  conditions,  roll  angles  will  not 
exceed  10  degrees  and  pitch  angles  are  not  expected  to  reach 
3  degrees.  Therefore,  the  following  assumptions  are 
justified  : 

sin  6  =0  (2-22) 

cos  0   =1  (2-23) 
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D.   CALCULATION  OF  VELOCITY  COMPONENTS 

The  various  foil  velocity  components  must  be  determined 
at  this  point  in  order  to  calculate  the  angles  of  attack  and 
side  slip.  The  total  velocity  components  of  each  foil  in 
body  axes  can  be  expressed  in  terms  of  the  craft  linear  and 
angular  velocities  by  considering  the  craft  geometry  as 
represented  in  Figure  7  and  Table  I. 
Center  Foil 


U   =  U  +  L    .  Q 
C         ZCF 


V   =  V  -  L    .P+L    .R 
C         ZCF       XCF 


N   =  W  -  L    .  Q 
C         XCF 


(2-24) 

(2-25) 
(2-26) 


Port  Foil 


U   =  U  +  L    .Q-L    .R 
P         ZPF       YPF 


V   =V-L    .P  +  L    .  R 
P         ZPF       XPF 


W   =  W  +  L    .P-L    .Q 

P         YPF       XPF 


(2-27) 
(2-28) 
(2-29) 


Starboard  Foil 


U   =U  +  L    .Q-L    .R 
S         ZSF       YSF 


V=V-L    .P  +  L    .R 
S         ZSF       XSF 


W=W  +  L    .P-L    .Q 
S         YSF       XSF 


(2-30) 
(2-31) 
(2-32) 


Mid-Foil  Segment 


0   =  U  +  L    .Q 
M         ZMF 


Y   =V-L    .P  +  L    .R 
M         ZMF       XMF 


(2-33) 
(2-34) 
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W   =  W  -  L    .Q 
H         XMF 


(2-35) 


•=Point  of  application  of  foil  forces 

0=Point   of   application   of   strut    forces 

L  =    L 

XCF  XCS 

L     =L     =  L     =L     =L 
XPF     XPS     XSF     XSS     XMF 

L     =  L     =  L     =  L 
YSF     YSS     YPF      YPS 

L     =  L     =  L 
ZPF     ZSF     ZMF 


L    and  L    are  variables  and  in  general  will  not  be 
ZPS      ZSS 


equal. 


c*  *cf 


FIGURE  7 
DEFINITION  OF  SYMBOLOGY 
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E.   ANGLE  OF  ATTACK  AND  ANGLE  OF  SIDE  SLIP 

The  angle  of  attack  of  a  foil  is  shown  in  Figure  8A  and 
is  defined  as  the  angle  whose  tangent  is  the  total  relative 
velocity  between  the  foil  and  water  in  the  bcdy  Z  direction, 
divided  by  the  total  relative  velocity  between  the  foil  and 
water  in  the  body  X  direction. 


a    =   arc  tan  (W   /0  ) 
r   r 


(2-36) 


If  the  foil  is  not  aligned  with  the  body  X  axis,  then 
the  foil  has  a  fixed  incidence  angle  and  the  total  angle  of 
attack  becomes  : 


a         =  a  +   a 

total      fixed 


(2-37) 


The  side  slip  angle  is  shown  in  Figure  8B  and  can  be 
represented  by  :  


/3   =  arc  tan 
i 


»  /  2     2    2 

V   /\/(U    +V   +W   ) 
ri  V   ri   ri   ri 


(2-38) 


During  actual  craft  operations  the  angles  will  be  less 
than  10  degrees,  so  that  the  inverse  trignometric  functions 
can  be  replaced  by  quotients.   Except  at  low  craft  speeds, 
the  forward  velocity  of  the  craft,  0,  is  much  greater  than 
the  velocities,  V  and  W.  Applying  these  simplifications,  the 
equations  become  : 

a   =  (tf   /U)  +  a 
i     ri       i-fixed 


0      -    V   /U 

i    ri 
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F.   EXPANSION  OF  FORCES  AND  MOMENT  TERMS 

With  the  aid  of  Figure  7  and  Table  I, the  forces  and 
moments  for  the  craft  are  easily  written.  It  is  important  to 
note  that  when  a  numerical  value  is  inserted  for  one  of  the 
alphabetical  dimensions,  the  sign  of  the  value  must  be 
positive  or  negative  according  to  positive  and  negative 
directions  from  the  axis  system  origin  at  the  C.G.  One  of 
the  assumptions  made  upon  commencing  simulation  was  that  the 
struts  contributed  with  zero  lift  and  that  the  foils 
contributed  with  zero  side  force.  Also  note  that  horizontal 
foils  make  no  contribution  to  the  body  Y  force  and  vertical 
struts  do  not  contribute  to  the  body  Z  force.   Therefore, 
all  body  Y  motion  originates  at  the  struts  and  all  body  Z 
motion  originates  at  the  foils. 

The  forces  and  moment  eguations  become  : 

F  =F    +F    +F     +F    +F    +F    +F     +mg  +T         (2-39) 
X   XCF   XPF   XSF    XMF  XCS       XPS  XSS  X   X 


F  =F    «-F    +F    +  mg  (2-40) 

Y   YCS   YPS   YSS     Y 

F  =F    +F    +F    +  F    +mg  (2-41) 

Z   ZCF   ZPF   ZSF    ZMF    Z 

L=  (F    L    )       +  (F  L   )   -(F  L  )   -  (F  L  )   -  (F  L  )  (2-42) 

Z  Y  PF    Z  Y   SF    Y  Z  PS    Y  Z  SS    Y  Z  CS 

M=-(F  L  )   -  (F  L   )   -  (F  L  )   -  (F   L  )   +  (F  L  ) 

Z  X  CF    Z   X  PF    Z  X  SF    Z   X  MF    X  Z  CF 

+  (F  L  )   +  (F  L  )   +  (F  L  )   +  (F  L   )   +  (F  L  ) 

X  Z  PF    X  Z  SF    X  Z  MF    X  Z   CS    X  Z  PS 

+  (F  L  )   +T  L  (2-43) 

X  Z  SS   X  ZT 

N=  (F  L  )   +  (F  L   )   +(F  L  )  (2-44) 

Y  X  CS    Y  X   PS    Y  X  SS 


All  of  these  forces  have  a  hydrodynamic  origin  except 
for  the  thrust  and  gravity  terms.  The  thrust  is  always 
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associated  with  the  body  axis  and  requires  no  further 
expansion.  The  gravity  terms  can  be  simplified  using  small 
angle  approximations  for  the  Euler  angles. 

The  standard  expression  for  hydrodynamic  forces  in 
water  axis  coordinates  is  : 

F   =  gAC       i=L,D,S  (2-45) 

i      i 


Where  L,D,S  represent  lift, drag  and  side  force  respectively. 
The  lift,  drag  and  side  force  coefficients  vary  as  functions 
of  angle  of  attack,  angle  of  side  slip,  submergence, 
velocity,  flap  deflections,  elevator  deflections  and  rudder 
deflections. 

Ideally,  mathematical  expressions  would  have  been 
developed  to  correctly  depict  the  interrelation  of  all 
variables  which  affect  the  force  coefficients.  However, 
except  for  the  conditions  when  relationships  are  linear, 
hydrofoil  technology  has  not  advanced  to  a  position  which 
would  yield  such  expressions.  The  difficulties  which  are 
encountered  in  deriving  a  mathematical  expression  include  : 

1)  The  ocurrence  of  cavitation  and  ventilation,  which 
are  not  completely  predictable. 

2)  The  even  more  unpredictable  cessation  of  cavitation 
and  ventilation. 

3)  The  nonlinearity  of  the  hydrodynamic  coefficients. 

4)  Lack  of  sufficient  test  data  to  completely  describe 
all  of  the  above. 
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Ill*   HYDROFOIL  SIMULATION 


A.   SIMULATION  OBJECTIVES 

The  objective  of  this  simulation  was  to  study  the 
dynamic  characteristics  of  motion  along  and  about  each  body 
axis.  That  is,  perturbations  should  be  introduced  to 
separately  excite  : 

1)  Motion  along  the  X-axis 

2)  Motion  about  the  X-axis 

3)  Motion  along  the  Y-axis 

4)  Motion  about  the  Y-axis 

5)  Motion  along  the  Z-axis 

6)  Motion  about  the  Z-axis 

The  following  assumptions  were  made  prior  to  beginning 
simulation  : 

1)  Craft  eguations  of  motion  are  valid  only  for  the 
foilborne  mode. 

2)  Weight  of  the  craft  remains  constant. 

3)  Hydrodynamic  coefficients  are  based  on  fully  wetted 
surfaces,  i.e.  no  cavitation  or  ventilation. 

4)  The  craft  may  be  operated  in  calm  water  or  in  a 
seaway. 
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B.   COMPUTES  SIMULATION 

The  computer  program  used  was  based   on   the   computer 
programs  of  fief . 1  and  Ref.2. 

As  in  the  preceding  References,  before  the  main  program 
could  be  assembled,  Fortran  subroutines  were  used  to  find 
the  hydrodynamic  coefficients  of  foil  lift  and  drag,  and 
strut  drag  and  side  slip.  The  SETUP  subroutine  stored  as 
data  points  the  simulation  curves  of  Ref.4.  The  subroutines 
INTERP  and  INTRP1  are  interpolation  subroutines  to  obtain 
the  proper  coefficient  from  the  above  simulation  curves.  For 
INTERP,  values  of  angle  of  attack  arid  forward  flap  angle 
were  used  to  enter  the  curves  and  obtain  a  value  of  foil 
drag.  Angle  of  side  slip  and  submergence  were  used  to  enter 
INTRP1  to  obtain  strut  drag.  Listings  of  the  subroutines  are 
shown  in  the  computer  program  section.  The  main  program  was 
then  assembled  using  Figure  9  as  a  reference  for  data  flow 
and  is  also  listed  in  the  computer  program  section.  In  one 
of  the  parts  of  this  main  program,  the  automatic  control 
block,  were  introduced  some  modifications,  needed  for  the 
main  objective  of  this  paper,  the  study  of  the  craft's 
turning  characteristics.  This  automatic  control  block  is 
still  open  for  future  studies. 
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C.   AUTOMATIC  CONTROL 

Control  signals  originate  at  the  pilots  controls,  at 
the  motion  sensors  and  at  the  position  transducer.  These  are 
combined  and  processed  in  the  electronics  to  produce  the 
signal  for  the  servo-valve,  The  pilot  will  be  able  to  move 
the  helm,  a  lever,  or  a  knob  to  introduce  steering,  altitude 
and  attitude  commands  ;  the  motion  sensors  sense  errors 
between  the  commanded  and  actual  craft  motions  and  produce 
electrical  signals  proportional  to  these  errors  ;  the 
position  transducer  provides  the  feedback  signal  so  that  the 
control  surface  stops  moving  when  it  has  reached  the 
position  corresponding  to  the  processed  error  signals. 

It  is  now  necessary  to  consider  how  the  control  system 
components  affect  the  equations  of  motion  through  the 
control  surface  deflections.  In  general,  a  signal  originated 
at  the  pilot  or  at  the  motion  sensors  will  be  modified  by 
the  individual  dynamics  of  the  components  in  the  signal  path 
before  the  control  surface  motion  actually  occurs.  Motion 
sensors  to  be  discussed  are  vertical  gyros,  rate  gyros, 
acceleromet ers  and  height  sensors. 
1  •  Vertical  G_yr o s 

A  vertical  gyro  here  will  produce  one  signal 
directly  proportional  to  the  pitch  angle  of  the  craft  and 
one  signal  directly  proportional  to  the  roll  angle  of  the 
craft. 

2-  Rate  G_yros 

A   rate  gyro  will  sense  an  angular  rate  about  one  of 
the  three  body  axes.  The  dynamics  of  the  gyro  are  such  as  to 
produce   an   output   signal   which   is  related  to  the  actual 
angular  rate  by  a  second  order  differential  eguation. 
3 .  Acceleromet ers 

Outputs  of  accelerometers  are  also  assumed  to  relate 
the  craft  motions  by  second  order  differential  equations. 

To  identify  the  accelerations  sensed  by  the 
instruments,  it  will  be  assumed  that  there  is  one 
accelerometer  at   the   C.G.,   mounted   such   as   to   measure 
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lateral  acceleration,  and  one  acceleroraeter  above  each  foil 
mounted  in  the  X-Y  plane  so  as  to  measure  accelerations 
parallel  to  the  Z-axis. 

The  lateral  acceleration  at  the  C.G.  is  : 

a     =  F  /m  (3-1) 

Y       Y 

Taking  angular  accelerations  into  account,  the 
accelerations  in  the  Z-direction  at  points  above  the  center, 
port  and  starboard  foil  are  : 

a    =  (F  /m)  -L    .Q  (3-2) 

ZC      Z      XCF 


a    =  (F  /m)  +L    .Q-L    .P  (3-3) 

ZP      Z      XPF     YPF 


a    =  (F  /ra)  +L    .Q+L    .p  (3-4) 

ZS      Z      XSF     YSF 

**•  Height  Sensor 

Ultrasonic  and  sonar  types  are  either  in  use  or 
expected  to  be  used.  For  this  study  it  will  be  assumed  that 
the  dynamic  characteristics  of  the  sensor,  whatever  the 
type,  will  be  such  as  to  give  a  direct  indication  of  the 
instantaneous  height  of  the  sensor  above  the  water. 
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3UDDER  ANGIE  -  -10* 


FIGURE      1 0 
SPAD^  RUDD1^   CHARAC'TRRISTICS 
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IV*  HYDROFOIL  TURNING  CHARACTERISTICS  STUDY 


A.   RUDDER  EFFECTS 


For  the  hydrofoil  steering  control,  the  struts  can  be 
viewed  as  vertical  toils  so  that  steering  can  be  achieved  by 
rotating  a  strut,  by  rotating  a  hinged  trailing  edge  or  by 
rotating  an  extension  of  a  strut  below  a  foil,  (spade  or 
ventral  rudder),  as  in  the  case  of  the  PC(H)-1. 

The  rudder  on  the  forward  strut  exerts  an  influence 
upon  the  forward  strut  side  force  coefficient.  Based  on 
Figure  10  the  necessary  changes  were  made  in  the  computer 
simulation  program.  To  document  the  rudder's  influence  in 
the  system  and  to  show  the  craft's  behavior,  several  figures 
are  presented,  corresponding  to  runs  using  different  speeds 
and  successive  rudder  angles  of  30,25,20,15,10,  and  5 
degrees  at  each  speed.  The  craft  started  moving  straight 
ahead,  the  rudder  angle  being  introduced  after  3  seconds, 
for  all  the  cases. 

1  •   Sp_eed=30  Knots 

For  a  forward  speed  U=30  knots,  Figure  11A  shows  the 
pitch  angle  of  the  craft,  which  had  an  oscillation  at  the 
beginning  of  the  run,  but  it  went  to  a  steady  state  value  of 
1.53  degrees,  with  the  craft's  bow  up,  in  4.5  seconds,  for 
all  the  six  different  rudder  angles  used.  The  result  was 
natural,  because  the  rudder  plays  a  negligible  role  in 
changing  the  value  of  the  foil  lift  coefficient  and 
therefore  in  the  craft's  pitch  motion. 

For  the  same  speed,  in  Figure  11B  the  different  roll 
angles  corresponding  to  different  rudder  angles  were 
represented.  When  the  ship  went  straight  ahead,  naturally 
the   roll  angle  was  0  degrees,  but  as  soon  as  a  ruddc      le 


of  30  degrees  was  introduced   after   3   seconds,   the   craft 

started   rolling,   going   to   a   steady  state  value  of  -1.29 

degrees  for  the  roll  angle  in  19.5  seconds.  On  the  other  end 

of   the   scale   for   a  rudder   angle  of  only  5  degrees,  the 

reaction  was   smaller,   the   initial   oscillation   was   also 

smaller  and  the  craft  went  to  a  steady  state  value  of  -0.216 

degrees  for  the  roll  angle  in   13   seconds.   Comparing   with 

Figure  4,  these  negative  values  show  a  roll  to  the  port  side 

direction,  although  the  ship  was  turning  to  starboard.   This 

result   would   be   perfectly   acceptable   for  a  conventional 

ship,  but  banking  in  the  opposite  direction  is   desired   for 

the  hydrofoil  because  of  the  high  speeds  of  operation. 

In  Figure  11C  were  represented  the   several   heading 

angles,   with  positive  values,  showing  that  the  craft  was  in 

fact  turning  to  starboard.  At   3   seconds   when   30   degrees 

rudder   angle   were   introduced,   the   craft  had  a  reaction, 

reaching  a   maximum   heading   angle   of   2.02   degrees,   but 

returning   to   a   steady  heading  value  of  1.39  degrees  in  23 

seconds.  For  a  rudder  angle  of  only  5   degrees   the   maximum 

heading   of   the   craft   was  0.33  degrees  ,  getting  a  steady 

heading  of  0.23  degrees  in  17.5  seconds,  because  the   rudder 

change  was  smaller  and  then  the  reaction  was  smaller. 

The   craft's   height   above   the    water    Z      was 

E1 

represented   in  Figure  11D.  Starting  with  a  value  of  0  feet, 

the  craft  went  to  3.9  feet  above  the  water  in  4.5    seconds. 

Comparing   with  the  dimensions  in  Figure  3A,  it  was  observed 

that  the  foils  were  well  inside   the   water.   Again   it   was 

observed  that  the  rudder  plays  a  negligible  role  on  the  lift 

motion  of  the  craft,  because  the  craft's   height   above   the 

water  was  the  same  for  all  the  rudder  angles. 

The  forward  flap  was  represented  in  Figure   11E.   It 

started   with   a  value  of  0  degrees,  reaching  a  peak  of  13.9 

degrees  to  help  the  transient  period.   After   5   seconds   it 

reached   a  steady  state  positive  value  of  10.8  degrees.  This 

positive  value  means  that  the  flap  was  oriented  in  the   down 
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direction  and  therefore  offering  more  resistance  to  the 
water,  provoking  more  lift  on  the  bow  part  of  the  craft.  It 
can  be  concluded  that  the  main  roles  of  the  forward  flap  are 
craft's  lifting  and  to  keep  a  stable  pitch  angle. 

In  Figure  11F  the  stern  mid  foil  segment  was 
represented.  After  a  transient  period  with  a  maximum  value 
of  -1.15  degrees,  it  reached  a  value  of  2.62  degrees  in  4.5 
seconds,  as  in  the  preceding  case  it  was  oriented  in  the 
down  direction  and  again  its  main  role  is  to  help  the  lift 
motion  of  the  craft. 

In  Figure  11G  the  stern  starboard  flap  was 
represented.  Its  value  was  0  degrees,  before  the  rudder  was 
introduced  after  3  seconds. Then  for  a  rudder  angle  of  30 
degrees,  the  flap  had  a  transient  period  with  a  maximum 
value  of  -2.33  degrees,  going  to  a  steady  state  value  of 
-1.86  degrees  in  16  seconds.  For  a  rudder  angle  of  only  5 
degrees, the  maximum  negative  value  was  only  of  -0.38 
degrees,  going  to  a  steady  state  value  of  -0.31  degrees  in 
13  seconds.  Having  the  starboard  flap  angles  different 
values  for  different  rudder  angles,  it  was  concluded  that 
the  stern  starboard  and  port  flaps  have  influence  in  the 
craft's  turns,  therefore  a  completely  different  conclusion, 
in  relation  with  the  conclusion  presented  for  the  forward 
flap  and  raid  aft  flap. 

On  the  next  groups  of  figures  the  same  elements  as 
before  will  be  represented.  Similar  conclusions  can  be 
applied,  therefore  only  the  data  obtained  will  be  reported. 
2.   Sp,eed=36  k n ot s ,  rudder  to  starboard 

The  pitch  angle  was  represented  in  Figure  12A.  It 
reached  a  steady  state  value  of  0.01  degrees  in  6.5  seconds. 

In  Figure  12B  were  represented  the  roll  angles  for 
the  different  rudder  angles,  showing  again  the  banking  in 
the  wrong  direction.  For  a  rudder  angle  of  30  degrees,  the 
roll  angle  had  a  maximum  negative  value  of  -1.29  degrees, 
going  to  a  steady  state  value  of  -1.12  degrees  in  15 
seconds.  For  a   rudder   angle   of   5   degrees,   the   maximum 


40 


negative   roll  angle  was  of  -0.21  degrees,  going  to  a  steady 
state  value  of  -0.18  degrees  in  12.5  seconds. 

The  heading  angles  represented  in  Figure  12C  reached 
higher  values  due  to  a  faster  speed,  increasing  gradually 
with  time,  meaning  that  with  the  increase  of  speed  the 
rudder  effect  on  the  craft's  turning  was  larger.  After  40 
seconds  of  run  time  with  30  degrees  of  rudder,  the  heading 
angle  reached  8.5  degrees  and  for  5  degrees  of  rudder,  the 
heading  angle  reached  1.42  degrees,  clearly  larger  values 
than  for  the  30  knots  speed. 

In  Figure  12D  the  height  above  the  water  was 
represented,  with  a  steady  state  value  of  1.65  feet  in  3.5 
seconds,  after  starting  on  0  feet. 

The  forward  flap  angle  was  represented  in  Figure 
12E.  On  the  transient  period  it  reached  a  maximum  value  of 
9.21  degrees,  going  to  a  steady  state  value  of  7.7  degrees 
in  5.5  seconds.  These  values  are  smaller  than  for  the  speed 
of  30  knots,  because  the  thrust  is  bigger  and  it  is  easier 
for  the  craft  to  lift. 

In  Figure  1 2F  the  stern  mid  foil  segment  was 
represented.  It  reached  a  maximum  value  of  -1.37  degrees  in 
the  transient  period,  jumping  to  a  steady  state  value  of 
1.74  degrees  in  7.5  seconds.  These  values  are  also  smaller 
than  for  the  30  knots  speed  by  the  same  reasons  of  the 
preceding  paragraph. 

In  Figure  12G  the  different  aft  starboard  flap 
angles  were  represented.  For  a  rudder  angle  of  30  degrees 
the  flap  angle  reached  a  maximum  value  of  -1.94  degrees, 
going  to  a  steady  state  value  of  -1.6  degrees  in  16  seconds. 
For  a  rudder  angle  of  5  degrees,  when  the  rudder  was 
introduced  the  maximum  value  was  -0.324  degrees,  reaching  a 
steady  state  value  of  -0.27  degrees  in  13  seconds. 
3 •   Speed=36  knots^  rudder  to  port 

With  the  same  speed  but  with  the  rudder  oriented  to 
the  opposite  board,  port,  there  were  no  changes  in  the  pitch 
mode,  the  only  changes  were  verified  with   the   roll   angle. 
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angle,  heading  angle  and  aft  starboard  flap  angle, 
represented  in  Figures  13A,  13B  and  13C.  The  corresponding 
figures  with  the  rudder  to  starboard  were  Figures  12B,  12C 
and  12G  and  comparing  these  two  groups  of  figures,  it  can  be 
observed  that  the  roll  angle,  heading  angle  and  aft 
starboard  flap  angle  have  exactly  the  same  values  but  with 
opposite  direction,  because  on  the  Y  body  axis  direction  the 
craft  is  perfectly  symmetrical. 
**•   S£eedf0.0  knots 

In  Figure  14A  the  several  roll  angles  for  different 
rudder  angles  were  represented. When  a  rudder  angle  of  3C 
degrees  was  applied  after  3  seconds  the  craft  started 
banking  in  the  wrong  direction  again,  having  a  transient 
oscillation  with  a  maximum  value  of  -1.18  degrees.  It  went 
to  a  steady  state  value  of  -1.032  degrees  in  15  seconds. 

In  Figure  14B  the  different  heading  angles  were 
represented.  The  craft  went  straight  ahead  and  started 
turning  when  the  rudder  angle  was  introduced  after  3 
seconds.  The  heading  angles  increased  gradually  with  time, 
for  a  rudder  angle  of  30  degrees  the  heading  angle  reached 
the  value  11.9  degrees  after  40  seconds  of  run  time.  For  a 
rudder  angle  of  5  degrees  the  maximum  heading  angle  reached 
1.99  degrees  in  40  seconds. 

In  Figure  1 4C  the  height  above  the  water  was 
represented.  It  went  to  a  steady  state  value  of  0.628  feet 
above  the  water  in  11  seconds. 

In  Figure  1 4D  the  forward  flap  was  represented. 
After  a  transient  oscillation,  with  a  maximum  value  of  6.9 
degrees,  it  dropped  to  a  steady  state  value  of  6.19  dgrees 
in  7  seconds. 

The  aft  starboard  flap  angles  were  represented  in 
Figure  14E.  For  30  degrees  rudder  angle,  the  flap  angle 
reached  a  maximum  negative  value  of  -1.78  degrees  and  went 
to  a  steady  state  value  of  -1.48  degrees  in  16  seconds.  For 
5  degrees  rudder  angle,  the  flap  angle  had  a  maximum 
negative  value  of  -0.29  degrees  and  went  to  a   steady   state 
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value  of  -0.25  degrees  in  11.5  seconds. 
5-   Sc-eedf 50  knots 

In  Figure  15A  the  different  roll  angles  were 
represented.  For  a  30  degrees  rudder  angle,  the  roll  angle 
had  a  maximum  value  of  -0.99  degrees  and  went  to  a  steady 
state  value  of  "0.88  degrees  in  10.5  seconds. For  a  5  degrees 
rudder  angle, the  roll  angle  reached  a  maximum  value  of 
-0.165  degrees  and  went  to  a  steady  state  value  of  -0.147 
degrees  in  8  seconds. 

The  heading  angles  were  represented  in  Figure  15B. 
For  a  30  degrees  rudder,  the  maximum  heading  angle  was  18.3 
degrees  in  40  seconds  of  run  time.  For  a  5  degrees  rudder, 
the  maximum  heading  angle  was  3  degrees  after  40  seconds. 

In  Figure  15C  the  forward  flap  angle  was 
represented.lt  went  to  a  steady  state  value  of  3.76  degrees 
in  3.5  seconds. 

In  Figure  15D  the  different  aft  flap  angles  were 
represented.  For  a  30  degrees  rudder  angle,  the  starboard 
flap  angle  reached  a  maximum  value  of  -1.51  degrees  and  went 
to  a  steady  state  value  of  -1.26  degrees  in  14.5  seconds. 
For  a  5  degrees  rudder  angle,  the  starboard  flap  angle 
reached  a  maximum  value  of  -0.25  degrees  and  went  to  a 
steady  state  value  of  -0.19  degrees  in  8  seconds. 
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B.   EFFECTS  OF  A  HELM  STEP  COMMAND 

Improving  and  completing  the  automatic  control  block  of 
References  1  and  2  and  based  on  the  block  diagram  of  Figure 
6.7  in  Ref.5,  a  simulated  automatic  control  system  was 
built,  including  the  dynamic  characteristics  of  the  sensors, 
the  dynamic  characteristics  of  the  compensation  filters, 
inputs  to  the  servo  valves  and  the  dynamic  characteristics 
of  the  hydraulic  actuators. 

There  was  an  important  limitation  and  perhaps  a  certain 
inaccuracy  in  the  values  used  as  corner  frequencies  of  the 
compensation  filters  and  the  values  used  as  gains  for  the 
control  actuation  servos,  because  References  1,2,3  and  4 
used  the  model  PC(H)->-1  for  their  studies  and  Ref.  5  used  a 
more  recent  model,  the  PKM.  Of  course  there  was  some  lack  of 
information,  but  some  interesting  results  and  some  data  were 
obtained  for  future  studies. 

As  mentioned  before, the  main  function  of  the  forward 
flap  was  to  lift  the  craft  and  keep  it  stable  above  the 
water  in  the  pitch-heave  mode.  To  confirm  this,  the  helm 
doesn't  actuate  the  forward  flap,  but  only  the  starboard  aft 
flap  and  the  port  aft  flap,  therefore  these  two  last  flaps 
will  be  responsible  for  the  craft's  turning  in  conjunction 
with  the  rudder  and  will  also  be  responsible  essentially  for 
the  craft's  rolling  and  finally  for  the  craft's  pitching. 
The  forward  flap  is  only  actuated  according  to  its  job  by 
the  forward  accelerometer  servo, by  the  height  sensor  servo 
and  by  the  pitch  angle  servo.  The  two  aft  flaps  are 
actuated  by  the  respective  accelerometer  sensors,  by  the 
pitch  angle  sensor,  by  the  roll  angle  sensor  and  by  the 
manual  helm  control.  The  rudder  servo  is  actuated  by  the 
roll  angle  servo,  by  a  yaw  rate  gyro  and  by  a  heading 
command  disengaged  or  engaged.  In  this  last  case  it  will 
stop  and  turn  back  the  rudder  when  the  desired  heading  is 
reached. 

Several   figures  will  be  shown  with   the   roll   angles, 
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heading   angles,   rudder   angles,   starboard  and  port  flaps 
angles  for  different  speeds  and  under  the  actuation  of   helm 
step   commands   of   150,120,90,60   and   30   degrees   to  the 
starboard  side,  at  each  speed. 
1 •  Speedf 30  knots 

In  Figure  16A  the  roll  angles  corresponding  to 
different  helm  angles  were  represented.  The  helm  was  moved 
to  starboard  after  12  seconds  of  straight  ahead  motion,  then 
the  craft  started  turning  in  the  right  direction,  also  to 
starboard.  For  a  helm  angle  of  150  degrees  the  roll  angle 
was  9  degrees  after  24  seconds  and  for  a  helm  angle  of  only 
30  degrees  the  roll  angle  was  3  degrees  after  30  seconds. 

In  Figure  16B  the  heading  angles  were  represented. 
The  craft  started  straight  ahead  and  when  a  starboard  helm 
was  introduced,  the  craft  started  turning  to  starboard, 
reaching  21.2  degrees  to  starboard  after  30  seconds  for  a 
helm  angle  of  150  degrees  and  only  3.4  degrees  to  starboard 
in  the  same  amount  of  time  for  a  helm  angle  of  30  degrees. 

Figure  16C  shows  the  different  rudder  angles  due  to 
a  helm  step  command.  While  the  ship  went  straight  ahead  the 
rudder  was  also  straight  ahead,  but  when  the  helm  was 
introduced  after  12  seconds, the  rudder  turned  to  starboard 
on  the  craft's  turning  direction.  Finally  when  the  helm  was 
removed  after  28  seconds,  the  rudder  went  back,  to  its  normal 
position.  For  a  helm  angle  of  150  degrees  the  rudder  reached 
a  value  of  5.1  degrees  to  starboard  and  for  a  helm  angle  of 
30  degrees  it  only  reached  0.78  degrees  to  starboard. 

In  Figure  16D  and  16E  the  port  aft  flap  and  the 
starboard  aft  flap  were  respectively  represented.  When  a 
helm  angle  was  introduced  after  12  seconds,  the  craft 
started  turning  and  banking.  After  those  12  seconds  the  port 
aft  flap  had  positive  values,  (down  direction)  and  the 
starboard  flap  had  negative  values,  (up  direction)  .  When  the 
craft  is  turning  to  starboard,  it  must  roll  to  starboard, 
therefore  a  lift  on  the  port  side  must  exist.  The  port  flap 
being  down,  when  the  craft  advances  more  force  is  exertad  on 
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that  flap  by  the  water  and  there  is  a  lift  force  on  that 
side.  The  starboard  flap  is  up  and  on  that  side  there  is  a 
down  force.  When  the  helm  goes  to  zero  degrees  after  28 
seconds,  the  craft  rolls  in  the  opposite  direction  and  it 
can  be  verified  by  the  figures  that  the  flaps  have  also 
moved  in  opposite  directions. 
2-   Speed=36  knots 

In  Figure  17A  the  roll  angles  corresponding  to 
different  helm  angles  were  represented.  When  the  helm  was 
introduced  after  10  seconds,  the  craft  started  rolling  to 
starboard.  For  a  helm  angle  of  150  degrees  the  roll  angle 
reached  6.2  degrees  in  18  seconds  and  for  a  helm  angle  of  39 
degrees  the  roll  angle  reached  2.55  degrees  after  27  seconds 
and  it  started  banking  in  the  opposite  direction  when  the 
helm  went  again  to  0  degrees. 

The  heading  angles  were  represented  in  Figure  17B. 
When  a  helm  angle  of  150  degrees  was  introduced  after  10 
seconds,  the  craft  started  turning  to  starboard,  reaching  a 
maximum  heading  angle  value  of  37.3  degrees  after  28 
seconds  and  for  a  helm  angle  of  30  degrees,  the  craft 
reached  a  maximum  heading  angle  value  of  6.9  degrees  after 
28  seconds. 

The  rudder  angles  were  represented  in  Figure  17C. 
For  a  helm  angle  of  150  degrees,  the  rudder  angle  reached  a 
value  of  5.37  degrees  to  starboard  after  24  seconds,  then  it 
started  going  back  in  the  opposite  direction  when  the  helm 
was  removed.  For  a  helm  angle  of  30  degrees  the  maximum 
rudder  angle  value  was  1.44  degrees  to  starboard  after  2'! 
seconds,  then,  as  before,  it  started  returning  in  the 
opposite  direction. 

In  Figure  17D  and  17E  the  port  and  starboard  flaps 
were  represented.  As  before  when  the  craft  was  banking  to 
starboard,  the  port  flap  had  positive  values  (down 
direction)  and  the  starboard  flap  negative  values  (up 
direction) .  Both  flaps  reversed  the  motions  when  the  craft 
rolled  back  to  its  normal  position  after  24  seconds. 
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3  •   Speed f  UP  Kl!0.!  s 

In  Figure  18A  the  roll  angles  for  different  helm 
angles  were  represented.  For  a  helm  angle  of  150  degrees, 
the  roll  angle  reached  a  value  of  9.38  degrees  after  19 
seconds.  For  a  helm  angle  of  30  degrees  the  roll  angle 
reached  a  maximum  value  of  2.43  degrees  after  23  seconds. 

In  Figure  18B  the  different  heading  angles  were 
represented.  For  a  helm  angle  of  150  degrees  the  heading 
angle  reached  26.3  degrees  in  24  seconds  and  for  a  helm 
angle  of  30  degrees  the  heading  angle  reached  12,12  degrees 
in  24  seconds. 

In  Figure  18C  the  different  rudder  angles  were 
represented.  For  a  helm  angle  of  150  degrees  the  rudder 
angle  had  a  maximum  value  of  5.83  degrees  to  starboard  after 
21  seconds  and  for  a  helm  angle  of  30  degrees,  the  rudder 
angle  had  a  maximum  value  of  1.16  degrees  to  starboard  after 
21  seconds. 

In  Figures  18D  and  18E  the  port  and  starboard   flaps 
were   shown.   Their   movements   were   the   same   as   before, 
according  to  the  craft's  rolling  motion. 
**•   Speed=50  knots 

In  Figure  19A  the  roll  angles  for  different  helm 
angles  were  represented.  For  a  helm  angle  of  150  degrees 
the  roll  angle  had  a  value  of  9.45  degrees  after  15  seconds 
and  for  a  helm  angle  of  30  degrees  the  roll  angle  had  a 
maximum  value  of  2.45  degrees  after  20  seconds,  then  it 
started  decreasing  because  the  craft  rolled  in  the  opposite 
direction  due  to  the  removal  of  the  helm  command  . 

In  Figure  19B  the  heading  angles  were  represented. 
For  a  helm  angle  of  150  degrees,  the  heading  angle  reached  a 
value  of  38.7  degrees  to  starboard  after  25  seconds.  For  a 
helm  angle  of  30  degrees  the  heading  angle  reached  a  value 
of  11.1  degrees  to  starboard  after  25  seconds. 

In  Figure  19C  the  rudder  angles  were  represented. For 
a  helm  angle  of  150  degrees,  the  rudder  angle  reached  a 
maximum  value  of  7.78  degrees  after  21  seconds  and   then   it 
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started  decreasing,  moving  in  the  opposite  direction  due  to 
the  removal  of  the  helm  command.  For  a  helm  angle  of  30 
degrees,  the  rudder  angle  reached  a  value  of  1.54  degrees 
after  21  seconds  and  again  it  started  decreasing. 

In  Figures  19D  and  19E,  the  port  and  starboard  flaps 
were  represented.  Their  motions  are  the  same  as  before,  but 
with  different  values. 
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C.   RESUL1S  WITH  THE  SPLITTING  OF  THE  FORWARD  FOIL 

The  forward  foil  was  split  in  two  parts  to  analyze  how 
it  improves  or  prejudices  the  craft1 s  pitching,  the  craft's 
rolling  and  the  craft's  turning  rate  in  the  desired 
directions. 

Some  changes  in  the  computer  simulaton  program  were 
necessary  to  achieve  this  trial.  Thus,  instead  of 
considering  a  unique  submergence  for  the  forward  foil,  two 
different  submergences,  Seep  and  Sees,  corresponding  to  the 
forward  port  foil  and  to  the  forward  starboard  foil  were 
considered.  These  two  foils  were  moved  independently,  thus, 
two  different  forward  flap  angles,  Elv1  and  Elv3  were 
considered.  Also  two  different  forward  foil  drag 
coefficients  and  two  different  forward  foil  lift 
coefficients  were  introduced.  Finally  the  block  diagram  of 
Ref.5  was  adapted  to  an  extra  forward  foil,  provoking  a 
small  change  on  the  simulated  automatic  control  system. 

With  a  helm  step  command  to  starboard,  the  craft 
started  turning  to  starboard  and  it  also  started  rolling  to 
starboard.  As  was  explained  in  part  B  of  this  chapter,  the 
two  forward  flaps  were  moved  in  a  manner  similar  to  that  of 
the  aft  port  flap  and  the  aft  starboard  flap,  thus,  the 
forward  port  flap  went  down  provoking  an  upward  lift  force 
on  the  port  side  and  the  forward  starboard  flap  went  up 
provoking  a  downward  force  on  the  starboard  side.  Several 
runs  with  different  speeds  and  different  helm  angles  were 
done.  In  a  general  way,  the  roll  angle  in  the  desired 
direction  was  improved  by  10%  and  the  turning  rate  was 
improved  by  8%,  but  the  pitch  angle  had  values  3%  smaller 
than  the  values  obtained  with  a  single  forward  foil,  because 
with  half  of  the  forward  foil  in  the  upward  direction,  there 
are  two  lift  forces,  (one  positive  and  the  other  negative) , 
prejudicing  the  craft's  bow  lift. 

Finally,  the  forward  flaps  were  moved  in  opposite 
directions  to  those  of  the  preceding   paragraph,   thus,   the 
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forward  port  flap  went  up  and  the  forward  starboard  flap 
went  down. The  results  were  worse  than  those  of  the  preceding 
paragraph. 
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D.   CALCULATION  OF  THE  CORRECT  ROLL  ANGLE 

Many  runs  with  the  craft  exposed  to  several 
disturbances  were  made.  Different  roll  angles  according  to 
the  existent  conditions  were  obtained,  but  now  a  question 
arises  :  are  these  roll  angles  sufficiently  near  the  optimum 
value  ? 

A  block  called  Calculation  of  the  Correct  Roll  Angle 
was  introduced  in  the  computer  simulation  program  to  answer 
this  guestion.  Knowing  the  craft's  turning  radius  for  a 
certain  speed  and  a  certain  helm  angle,  the  craft's 
centrifugal  force  can  be  calculated.  The  craft's 
gravitational  force  in  the  Z  direction  was  already  known, 
then,  the  force  perpendicular  to  the  deck,  called  net  force, 
can  be  calculated,  by  the  gravitational  and  centrifugal 
forces.  The  angle  between  the  gravitational  force  and  the 
net  force  perpendicular  to  the  deck  is  the  correct  roll 
angle. 

In  Figures  20A  and  20B  the  craft's  trajectories  in  the 
X-Y  plane  for  a  speed  of  40  knots  and  helm  step  commands  of 
20  and  60  degrees  respectively  were  represented.  The  craft's 
turning  radius  was  calculated  at  the  point  where  the  beading 
angle  was  reversed  180  degrees  in  relation  with  the  initial 
heading  angle  at  the  beginning  of  the  run,  since  this  will 
be  approximately  the  steady  state  value. 

In  Figure  20C  the  craft's  turning  path,  for  a  helm 
command  of  160  degrees,  a  speed  of  50  knots  (extreme 
conditions)  and  a  run  time  of  156  seconds  was  represented. 
The  turning  path  was  approximately  a  circle,  with  an 
approximate  turning  radius  of  410  feet. 

The  values  of  the  actual  roll  angle  matched  closely 
with  the  values  of  the  correct  roll  angle,  for  the  steady 
state  period,  although  there  were  differences  during  the 
transient  period.  These  values  are  represented  in  Table  II. 

In  Figure  20D  the  craft's  turning  path  for  a  helm  angle 
of  160  degrees  and  a  speed  of  30  knots  was  represented.   The 
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turning  radius  was  smaller,  approximately  300  feet. 

The  values  of  the  actual  roll  angle  matched  again  with 
the  values  of  the  correct  roll  angle,  for  the  steady  state 
period.  These  values  are  shown  in  Table  III. 


TABLE  II 


Roll  angles  for  helm=160  degrees  and  speed=50  knots 


Time 

Actual  roll  angle 

Correct  roll  angle 

(Seconds) 

(Degrees) 

(Degrees) 

17.3 

13.98 

5.65 

26.0 

17.08 

10.22 

34.7 

18.33 

14.07 

43.  4 

18.86 

18.  19 

52.0 

19.14 

19.31 

60.7 

19.27 

19.39 

70.3 

19.36 

19.46 

78.  1 

19.39 

19.52 

86.8 

19.39 

19.53 

95.4 

19.39 

19.53 
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TABLE  III 


Roll  angles  for  helm=160  degrees  and  speed=30  knots 


Time 

Act 

ual  roll 

an 

.gle 

Correct  roll  angle 

Seconds) 

(Degrees 

') 

(Degrees) 

19.  2 

10.77 

7.32 

28.9 

14.53 

10.51 

38.5 

14.88 

13.84 

48.2 

15.29 

14.59 

57.8 

15.47 

15.12 

66.9 

15.47 

15.72 

75.2 

15.47 

15.72 
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V«   CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  STUDIES 


The  system  was  made  stable  in  six  degrees  of  freedom, 
with  the  following  limitation  :  the  dimensions,  the  foil 
areas,  the  weight  and  the  thrust  of  the  model  PC(H)-1  were 
used  in  References  1  and  2  and  in  the  computer  simulation 
program.  Due  to  a  lack  of  information,  because  only  very 
general  information  was  available,  certain  values  of  the 
automatic  control  system  of  the  model  PHM,  instead  of  the 
PC(H)-1,  had  to  be  used. 

A  study  of  the  craft  turning  characteristics  was  made 
in  this  thesis.  The  simulated  automatic  control  system  was 
improved  and  adapted  to  perform  this  study. 

Roll  was  controlled  by  differential  displacements  of 
the  aft  flaps,  actuated  by  the  helm  control,  while  pitch  was 
essentially  controlled  by  the  forward  flap. 

The  rudder  has  negative  effects  on  the  hydrofoil's 
rolling  direction,  which  should  be  the  same  as  the  turning 
direction,  due  to  the  high  speeds  of  operation. 

It  was  not  worth  while  to  split  the  forward  foil. 

Recommendations  for  further  studies  include  : 

1.   Explore   more   deeply    the    advantages    or    the 


disadvantages 


the   splitting   of   the   forward   foil, 


combining  it  with  the  motions   of   the   aft   flaps   and   the 
rudder. 

2.  Study  the  response  of  the  craft  in  a  sinusoidal  sea 
and  in  a  random  sea.  Design  improved  control  systems  for  sea 
state  operation. 

3.  Simulate  the  transformation  from  hullborne  operation 
to  foilborne  operation. 
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APPENDIX  A 


SUMMARY  OF  COMPLETE  EQUATIONS 


EQUATIONS  FOR  ACCELERATION  IN  BODY  AXIS 

U  =  |fx+RV-QW  (A-l) 

V  =  -F  +PW-RU  •                (A- 2) 

W  =  -F^+QU-PV  (A-3) 
in  /ii 

P  =  T^-[L-QR(Izz-IyY)+(R+QP)Ixz]  (A-4) 

Aa 

r  =  I^-[N-PQ(Iyy-Ixx)-(QR-P)Ixz]  (A-5) 


VELOCITY    IN   BODY   AXES 

U   =   /   Udt  (A-6) 

V  =   /   Vdt  (A-7) 

W  =   /   Wdt  (A-8) 

P  =   /   Pdt  (A-9) 

"     Q  =   /   Qdt  (A-10) 

R  =   /    Rdt  (A-ll) 

TRANSFORMATION   OF   VELOCITIES    FROM   BODY   TO  EARTH  .AXES 

U-n,=Ucos0cosijj+V(cos!];sinesin4)-sin^cosc{)  )+Y/(costjjsinQcos({)  +  sini;'Sin£ ) 

*  (A-12) 

V   =UcosGsiniJ;+V(cosij;cos4>  +  sinijjsin9sin<J>  )+W(sini|jsin9cos(j)-cos^sin<)) 

*  (A-13) 

W   =-UsinQ+VcosOsin<{>+WcosGcos<{>  (A- 14) 
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$   -  P+t^sinG 

0    =  Qcos<J>-Rsin<{) 

ty   =    (Qsin<J>+Rcos<{>)cosO  +  (|-P)sin0 

POSITION  IN  EARTH  AXES 


XE 

— 

/  U£dt 

YE 

= 

/  vEdt 

.  ZE 

= 

/  W„dt 
h 

* 

= 

f   Jdt 

e 

= 

/  6dt 

* 

= 

/  }dt 

EXPANSION  OF  FO 

Fx  = 

FY  = 

F   - 
Z 

L  = 

M   = 
N   = 


;a-i5) 
;a-i6) 

A-17) 


E  F    +  E  F    +  ms  +T 

L  XIF    L      XIS     ^X   X 

Z  FYIS  +  mgY 

1  FZIF  +  mgZ 

(FZLY)PF+(FZLY)SF"(FYLZ)PS"(FYLZ)SS~(FYLZ)CS 

E  (FZLX)IF+Z  (FXLZ)IF+E  (FXLZ)IS+TXLZT 

E  (FyLx)IS 


EXPANSION  OF  GRAVITY  TERMS 


mgx  =  -mg  sine 


mgy  =   m5   cos6sin(|) 
mg      =  mg   cosGsiiK}) 


:a-i8) 

;a-i9) 

A-20) 
A-21) 
A-22) 
;A-23) 


:A-24) 
A-25) 
:A-26) 
A-27) 
;A-28) 
;A-29) 


;a-30) 
;a-31) 
;a-32) 
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FOIL  VELOCITY  COMPONENTS  IN  BODY  AXES 
Center  Foil: 

Uc  =  U+LZCF-Q  (A-33) 

vc  =  v-lzcf'p+lxcf*q  <a-34> 

Wc  =  W-Lxcp-Q  (A-35) 

Port  Foil: 

up  =  ^zpf'^^pf'11  (A"36) 

VP  =  V-LZPF*P+LXPF*R  (A~37) 

Wp  =  W+LypF-P-LxpF-Q  (A- 38) 

Starboard  Foil: 

us  =  u+lzsf'q"lysf"r  (A_39) 

VS  =  V-L2SF-P+LXSF-R  (A"40) 

WS  =  W+LYSF-P-LXSF-Q  (A'41) 

Mid  Foil: 

UM  =  U+LXMF,Q  <A"42) 

VM  =  V-LZUF-P+LXMF-R  <A~43> 

WM  =  W"LXMF-Q  <A"44> 
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TRANSFORMATION  OF  WATER  PARTICLE  ORBITAL  VELOCITY  FROM  EARTH 
TO  BODY  AXES 

UWI  =  UE^icos^cos0+VE\visin^cose~wE^isine  (A-45) 

VWT  =  UFWT  (cos^sin9sin(J)-sirn|;cos({)  )+Vv,TrT  (sin^sin9sin(J)+cosi|;cos(J)  ) 

+WEWIcos6siri({)  (A-46) 

WWI  =  U  WI(cosi|;sin9cos({)+sin^sin(J)  )+VWFI(sin<j;sin0cos(j)-cos^sin<J> ) 

+WEWIcos9cos<J>  (A-47) 

Each  of  these  equations  must  be  repeated  for  each  foil. 

RELATIVE  VELOCITY  COMPONENTS 

URI  =  Uj  -  UWI  (A-48) 

'  VRI  =  VI  "  VWI  (A"49> 

WRI  =  Wj  -  WWI  (A- 50) 

One  set  of  equations  for  each  foil. 

ANGLES  OF  ATTACK  AND  SIDE  SLIP 

•   QI  =  arCtan  WRl/URI  +  aI-FIXED  (A~51)* 

VRT 

BT  =  arcsin  RI    —  (A-52)* 

/U2  +V2  +W2 
'URI   RI   RI 

*  One  set  of  equations  for  each  foil.   Angles  of  attack  and 
side  slip  are  calculated  in  radians. 
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TOTAL  RELATIVE  VELOCITY  AT  A  PARTICULAR  FOIL  OR  STRUT 


V  =  -%+V^+W^  (A-53) 


HYDRODYNAMIC  FORCES  IN  WATER  AXES 

(Each  foil  equation  adn  strut  equation  must  be  repeated  for 

foil  and  strut,  respectively) 

Foil  Lift: 

FL  =  £pV2AFCF  (A- 54) 

Foil  Drag: 

FDF  =  -PV2AFCF  (A"55) 
Strut  Drag: 

FDS  "  *PV>ASCDS  <A"56> 
Strut  Side  Force: 

Fs   =  ipV2AgCs  (A-57) 

TRANSFORMATION  OF  FOIL  FORCES  FROM  WATER  AXES  TO  BODY  AXES 

FyTF  =  -FTDFcosacosB+F- TFsina  (A-58) 

FYIF  =  -FDIFSinB  "   (A"59) 

FZIF  =  -FDIFsinacosB-FLIFcosa  (A-60) 
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//    EXEC    DSLtPARM.C='NOSOURCE,NOMAP' 

//DSL.  INPUT    DD    * 

D  CCMMON    MATRIX(539) 

CONTRL     FINTIM=3O.O,DELT=0.025,DELS=0.1 

INTEG    RKSFX 

INTGER     IFLAG, IFLAGl,NPLOT, KUM,NCUR 

CONST    NPL0T=1  ,NCUR=1, IFLAG=0 

CONST    WEIGHT=262527,     IXX=601200,     I YY  =  '+286000, . . . 

IZZ=4l 740QO,G=32.2 

CONST  ACF=3  2.7,APF=36.3,ASF=36.8,  AMF=76.0,... 

PI=3. 141593 

CONST    LXCF=41 .5  9?LXCS  =  41 . 59, LXPF=-1 6 . 9  8 , LXPS=-  16 . 98, . . . 

LXSF=-16.93,LXSS=-16. 9  8 , LXMF=-16 . 98, L YCP=-5. 0 , LYCS  =  5.0, . .  . 

LYPF=-12.42,LYPS=-8.0,LYSF=12.42,LYSS=8.O,LZCF=15.5,... 

LZCS=13.25,LZPF=17.14,LZPS=14.0,LZSF=L7.14,LZSS=14.0,... 

LZMF=17.14,LXT=-16.9,LZT=17. 7,LXH=6  8.0 

CONST    Kl=5.0,K2=-96.0,K3=-82.0 

CONST    K4=1.5,K5=10.0,K&=3.5,K7=1.42,K8=l.5,K9=-2.0,... 

K10=0.44,K11=-13.0,K12=3.0,K13=3.0,K14=0.014 

CONST    RST=0.0 

STORAG     ICK3)  ,  NU1  (  3  )  ,  DEN1  (  4)  ,  IC2  (  2)  ,  NU2(  2)  ,DEN  2(  3  ) ,  .  .  . 

IC3(3) ,NU3(2) ,DEN3<4) 

TABLE    IC1(1-3)=3*0.0,NU1( 1-3 )  =  1 .0  ,3. 0 1 ,0.03, . .  . 

DEN  I (1-4) =1.0, 110. 1,3011.0,300.0, I C2(  1-2)  =  2-0.0,. .. 

NU2(1-2)=1.0,1.0,CEN2( 1 -3 ) = 1 .0 , 8  7 . 7, 1 8  85 .0 , IC 3 ( 1-3) =3*0 . 0, . . 

NU3( 1-2) =30. 0,5 8.5, DEN3( 1-4) =1.0, 30. 25, 0.0 15, 0.45 

PARAM    UZERO=40.0 

PARAM    RUDR=O.C 

PARAM    HELMR=2.79 

PREPAR     .  1,RTIME,PHI1,PSI  1,XE,YE 

GRAPH    XE,YE 

PRPLCT 

LABEL    RESPONSES    WITH    HELM    STEP    COMMAND 

PRINT    0. 5, R TIME, PHI ,PHI3,PHI1 ,PHI4, PS  I , PS I 1,XE,YE ,UE,VE, . .. 

RGIR,FGRAZ,FCENT,FNET 

INITIAL 

CENSE  =  2  .0 

UINIT=UZER0*2000. 0*3. 0*1. 0/3600.0 

U=UINIT 

TX=3960.0 

ALM=0. 074269 

ALS=0. 074269 

ALP=0. 074269 

ALF=0. 037111 

PREF=0.0 

TRFF=0.0 

THETA2=0.0 

PHI2=0.0 

PSI2=0.0 

MASS=WEIGHT/G 

CCNV=PI/180.0 

HREF=11  .5 

HREF2  =  -11  .5 

LENGTH=115.75 

FLAG=0.0 

IFLAG1=0 

EPS=25.0 

IF(IFLAG.EQ.O)    CALL    SETUP(FLAG) 

IFLAG=1 

IF(FLAG.EQ.l)     CALL    ENDRUN 

DERIVATIVE 

NO SORT 

*      ANGLES  OF  ATTACK  AND  SIDE  SLIP 

VRC=V-LZCF*P+LXCF*R 
VRS=V-LZSF*PU.XSF*R 
VRM=V-LZMF*P+LXMF*R 
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SORT 


VRP=V-LZPF*P+LXPF*R 
WRP=W+LYPF*P-LXPF*Q 
KRS*W+LYSF*P-LXSF*Q 

WFM=W-LXMF*Q 
WRC=W-LXCF*Q 
ALFC=WRC/U+ALF 
ALFAP=WPP/U+ALP 
ALFAS=WRS/U+ALS 
ALFAM=WRM/U*ALM 
BETAC=VRC/U 
BETAP=VRP/U 
BETAS=VRS/U 
* 

*  CALCULATION    OF    EULER    ANGLES 

PHID1=P 

TFETD1=G 

PSID1=R 

THETA=INTGPL(THETA2,THETDi) 

PHI=INTGRL(PHI2,PHID1) 

PSI  =  INTGRL(PSI2,PSI01  ) 

WE=-U*SIN{ THETA) +V*COS(THETA)*SIN(PHI)+rt*COS(THETA) 

*COS(PHI) 

ZE=INTGRL(HREF2,WE) 

ZE1=ZE-HREF2 

NOSORT 

*  SUBMERGENCE    TERMS 

SECF=ZE-LXCF*THETA+LZCF*COS(PHI  ) 

SEPF=ZE-LXPF*THETA+LYPF*SIN{PHI ) +LZPF*COS ( PHI ) 
SESF=ZE-LXSF*THETA*LYSF*SIN( PHI ) +LZSF*COS ( PHI ) 
SEMF=ZE-LXMF*THETA+LZMF*COS(PHI) 

*  HYDRODYNAMIC  COEFFICIENTS 

HELMR1=STEP(RST) 

HELM=HELMR*HELMR1 

RUD11  =  K7*PHH-K9*R 

RUD=REALPL(0.0,0.159,RUD11) 

ELVl=ELV*CONV 

ELVD=LIMIT(-5.0t5.0,ELVl) 

FLAD=FLAP*CONV 

FDEG=LI MIT(-5.0,5.0,FLAD) 

ARPD=ARP*CONV 

ARP1=LIMIT(-5.0,5.0,ARPD) 

ARSD=ARS*CONV 

ARS11  =  LIMIT<-5.0,5.0,AP.SD) 

CALL     INTERPC ALFC,ELVD,COEFC,2) 

CALL     INTERP(ALFAP,ARP1  ,COEFPX,3) 

CALL     INTERP( ALFAS,ARS11 ,C0EFSX,3) 

CALL     INTERP( ALFAMt FLAP,COEFMX, 1  ) 

CALL     INTRP1<BETAC,SECF,CFDS,2,THETA,THETA1 ) 

CALL     INTRPK  BET AP , SEPF ,CPDS,  I  ,  THETA ,  THE  T  Al } 

CALL     INTRPK  BET  AS, SESF ,CSD Sf 1 1 THETA , THETA 1 ) 

*  FORCES    IN    THE    WATER    AXES 
* 

FGRAX=-MASS*G*S IN (THETA) 
FGRAY=MASS*G*SIN(PHI) 
FGRAZ=MASS*G*COS<PriI i 
CCS=-0.043  75*BETACK).04375*RUD 
ACCS=CCS*4.25*SECF 
FSO0.5*DENSE*<U**2l*ACCS 
FDSC=0.5*OENSE*( U**2) *CFDS 
PCS=-2. L5*BETAP 
APCS=PCS*4.25*SEPF 
FPS*'0.5*DENSE*(U**2J*APCS 
FDPS=0. 5*DENSE*( U**2) *CPDS 
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SCS=-2.15*8ETAS 
ASCS=SCS*4.25*SESF 
FSS=0.5*DENSE*( U**2)-ASCS 
FCSS=0.5*DENSE*(U**2) *CSDS 
CLCF=5. l5*ALFC+0  .036*ELV  +  0. 13 
FLC  =  0.5*DEi\lSE*(U**2l *ACF*CLCF 
FDC=0.5*DENSE*(  U**2)  *ACF*COEFC 
CLPF=4.iO*ALFAP+0.018*ARP 
FLP=0 . 5  *  DENS  E*  <  U#*2 ) * A  P  F*C  L  PF 
FDP=0.5*DENSE*(U**2)*APF*COEFPX 
CLSF=4. 10*ALFAS+0.018*ARS 
FLS=0.5*DENSE*(U**2)*ASF*CLSF 
FDS=0.5*DENSE*(U**2)*ASF*COEFSX 
CLMF=5. 15*ALFAM+0.036*FLAP 
FLM=0.5*DENSE*(U**2)*AMF'*CLMF 
FDM=0.5*DENSE*< U**2)*AMF*COEFMX 
FLCZ=0. 5*DENSE*(U**2)*ACF*CLCF 
FCCZ=0. 5*DENSE*(U**2) *ACF*COEFC 
FLPZ=0. 5*DENSc*  (  U**2)  *APF*CLPF 
FCPZ=0.5*D£NSE*(U**2)*APF*COEFPX 
FLSZ=0. 5*DENSE*< U**2) *ASF*CLSF 
FDSZ=0. 5*DENSE*< U**2) *ASF*COEFSX 
FLMZ=0.5*DENSE*(U#*2) *AMF*CLMF 
FDMZ=0.5*DENSE*(U**2) *AMF*COEFMX 

TRANSFORMATION    TO    BODY    AXES 

FXCF=-FDOALFC*FLC 
FXPF--*FDP+ALFAP*FLP 
FXSF=-FDS*ALFAS*FLS 
F  X  M F=  -F  DM  +  AL  F  A  M *  F  L  M 
FXCS=-FDSC-BETAC*FSC 
FXPS=-FDPS-BETAP*FPS 
FXSS=-FCSS-BETAS*FSS 
FYCS=-BETAC*FDSC+FSC 
FYPS=-BETAP*FDPS+FPS 
,  FYSS  =  -BETAS*FDSS<-FSS 

FZCF--FDC*ALFC-FLCZ 
FZPF=-FDPZ*ALFAP-FLPZ 
FZSF=-FDSZ*ALFAS-FLSZ 
F  Z  MF  =-F  DM  Z  *  A  L  F A M- F L  M  Z 
LZCS=LZCF-0.4*SECF 
LZPS=LZPF-0.4*SEPF 
LZSS=LZSF-0.4*SESF 

FXCNE=FXCF+FXPF+FXSF*FXMF+FXCS*FXPS+FXSS 

*  SUMMATION    OF    FORCES    AND    MOMENTS 

FX=FXONE+FGRAX+TX 

UDOT=(l ./MASS)*FX-Q*W 

IFUFLAG1.EQ.O)     UDOT=--0.0 

FY=FVCS+FYPS+FYSS+FGRAY 

FZ=FZCF+FZPF+FZSF+FZMF+FGRAZ 

l  =  FZPF*LYPF  +  FZSt=*LYSF-FYPS*LZPS-FYSS*LZSS-FYCS*LZCS 

MFWD=FZCF*LXCF 

N=-MFWD-FZPF*LXPF-FZSF*LXSF-FZMF*LXMF#-.  .  . 

FXCF*LZCF+FXPF*LZPF+FXSF*LZSF+FXMF*LZMF* FXCS... 

*LZCS*FXPS*LZPS+FXSS*LZSS+TX*LZT 

N=FYCS*LXCS+FYPS*LXPS+FYSS*LXSS 

*  VELOCITY    INTEGRALS 
# 

SORT 

■if. 

VDCT=( l./MASS)*FY-R*U 
WDOT=( i ./MASS)*FZ+Q*U 
PDCT=(1 .0/IXX)*L 
CDOT=(  I  .0/IYY)~  '1 
RDCT=( 1  .0/IZZ)*N 
U=INTGRL(UINIT,UDOT) 
V=INTGRL(0.0,VDOT) 


IIP 


W=INTGRL(0.0,WDOT) 

P=INTGRL(0.0,POOT) 
Q=INTGRL(0.0,QDOT) 
R=INTGRL(0.0,RDOT) 
# 

*  SENSED  ACCELERATIONS 

AZC=FZ/MASS-LXCF*QDOT 

AZS=FZ/MASS+LXSF#QOOT+LYSF*PDOT 

AZP=FZ/MASS+LXPF*QDOT-LYPF*PDOT 

*  AUTOMATIC  CONTROL  BLOCK 

*  * 

HACT=-ZE+LXH*THETA 

DELH=HREF-HACT 

DELHC=REALPL( 0.0,0. 04, DELH) 

DETH=TREF-THETA 

FLAP=K2*DETH 

B=AZCfK14*(PHl**2) 

AZCC=TRNFR(2,3,IC3,NU3,DEN3,BI 

AZSC=R£ALPL( 0.0,0. 033, AZS) 

AZPC=REALPL( 0.0,0.033, AZP) 

THETAC  =  TRNFR(2f 3 , I C 1 , NU1 , DEN1 , 1000 .0*TH ETA ) 

PHIC=TRNFR(1 ,2,IC2,NU2,DEN2,1865„0*PHI  ) 

HELMC=REALPL(0.0, 1.0, HELM) 

ELV11=K11*THETAC+K12*THETA+K13*DELHC+K4*AZCC 

ELV=REALPL(0.G,0.159,ELV11) 

ARS2=K4*AZSC*-K5*THETAC+K6*PHIC-K10*HELMC 

ARP2=K4*AZPC+K5*THETAC-K6*PHIC*K10*HELrtC 

AFS=  REAL  PL  (CO,  0.1 59, ARS2) 

ARP=REALPL(0.0,0.15  9,ARP2) 

ThETA3=(THETA-THETA2)*180.0/PI 

PHI  1= (  PHI-PHI  2) -180. O/PI 

PSI1=< PSI-PSI2)*180.0/PI 

RUD1-RUD*I80.0/PI 

*  CALCULATION  OF  THE  CORRECT  ROLL  ANGLE 

UE=U*CCS(THETA)*COS(PSI)*V*(COS(  PS  I J  *S  IN  (  THETA)  *.  .  . 

SIN(PHI  J -SIN (PS  I  )*COS(PHI  )) +W*( COS (PS  II *SI N ( THETA ) .. 

*CCS(PHI )+SIN(PSI )*SIN(PHI  )  ) 

VE  =  U* CO S< THETA) *SIN(PSI )*V*(COS( PSI) -COS ( PH I ) ♦  .  .  „ 

SIN(PSI ) *S IN ( THETA )#S IN(PHI  ) ) +W* ( S IN ( P S I )*S  IN (THETA) 

*CCS( PHI ) -COS (PS  I  )*SIN(PHI ) ) 

YE=INTGPL(0.1,UE) 

XE=INTGRL<0.1 ,VE) 

RGIR=(  (  XE**2)<-<  YE** 2)  )/XE 

FCENT=MASS*( ( ( UE**2) +( VE**2 ) l/RGIR) 

FNET=SQRT  ( (FGRAZ**2  )  f  (  FCE.MT**2  )  ) 

PHI3=ATAN2( FCENT, FGRAZ ) 

PHI4=PHI3*18G.0/P I 

*  CALCULATION  OF  REAL  TIME 
PTIME=( LENGTH/UINIT)*TIME 

DYNAMIC 

SAMPLE 

CALL  DRWGdt  ltXEtYE) 
TERMINAL 

CALL  END RW(N PLOT) 
END 
STOP 
FORTRAN 
C 
C 

C  SUBROUTINE  SETUP 

C      THIS  SUBROUTINE  FORMS  THE  DATA  POINTS  COMPRISING 
C      THE  CURVES  OF  FOIL  AND  STRUT  DRAG  COEFFICIENTS 
C      AND  SIDE  SLIP  COEFFICIENTS  INTO  MATRACES. 
C 
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2 

100 


10 


11 

12 


SU 
CC 
DI 
1CF 
RE 
IC 
IF 
IF 
MA 
J  = 
K  = 
IC 
RE 
F0 
II 
IF 
GC 
CC 
DO 
J  = 
IF 
IF 
CC 
IC 
IC 
IF 
CC 
GO 
II 
IF 
J  = 
K  = 
IC 
MA 
IF 
DO 
J  = 
IF 
IF 
CO 
IC 
IC 
IF 
CC 
GC 
I  1 
IF 
J  = 
K= 
IC 
MA 
IF 
DO 
J  = 
IF 
IF 
CC 
IC 
IC 
IF 
CC 
GO 
I  1 
IF 
J  = 
K  = 
IC 
MA 
IF 
DC 
J  = 


GO  TO  51 
8 ) , MATRIX 


BROUTINE  SETUP(FLAG) 

MMCN  COEFMXtCOtFCX,COEFOX,CFDSl  ,CADS1 

MEN'S  I  ON  C0EFCX(3T33)  ,COEFOX(  i  T  33)  ,  COEF  MX  (  3 

DS1 ( 11,11)  , CADS  if  11  ,11) , INF  0(8) 

AL  INFO 

=  0 

LAG  =  0 

LAG1=0 

TRIX=1 

0 

1 

CUNT=0 

AD(5,100)  INFO 

RMAT(7F10.8) 

=  1 

(IFLAG1.EQ.1) 

TO  (4,5,6,7, 
NTINUE 

1  1=1,8 
J  +  l 

(J.GT.3)  K  =  K  +  1 
(J.GT.3J  J=l 
EFOX( J5K)=INFO(  I ) 
CUNT=IC0UNT+1 
=  IC+1 

(IC0UNT.EQ.99)  GO 
NTINUE 

TO  2 
=  1 

(II. LE 
0 
1 

CUNT=0 
TR1X=2 
(I1.GT.8)  GO  TO  2 

9  1=11,8 
J  +  l 

(J.GT.3)    K=K+1 
(J.GT.3)    J=l 
EFCXC J,KJ  =  INFO{ I  ) 
0UNT=ICQUNT+1 
=  IC  +  1 


33) 


TO    3 


)    11=11+1 


(  IC0UNT.EQ.99)     GO    TO    10 
NTINUE 

TO    2 
=  1 

(I1.LE.8)     11=11+1 
0 
1 

CUNT=0 
TRIX=3 
( I1.GT.8)  GO  TO  2 

11  1=11,8 
J+l 

(J.GT.3I  K=K+1 
(J.GT.3)  J=i 
EFMX( J,K)=INFO( I ) 
CUNT=IC0UNT+1 
=  IC  +  1 

( IC0UNT.EQ.99)  GO  TO  12 
N'T  I  NU  E 

TO  2 
=  1 

(I1.LE.8)  11=11+1 
0 
1 

CUNT=0 
TRIX=4 
(I1.GT.8J  GO  TO  2 

13  1=11,3 
J  +  l 
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8)     11=11+1 
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FORTR 
C 

c 

c 
c 


IF(J.GT.ll)    K  =  K  +  1 

IFU.GT.ll)     J=l 

CADSK  J  ,KI=JNF0<  I) 

ICCUNT=ICGUNH-1 

IC=IC+1 

IF( IC0UNT.EQ.121)  GO 
13  CONTINUE 

GC  TO  2 
15  11=1 

IF( II. LE 

J  =  0 

K=l 

ICCUNT=C 

MATRIX=5 
8   I  F ( 1 1 .GT.8)  GO  TO  2 

DC  14  1=11,8 

J  =  J  +  1 

IFU.GT.ll)  K=K  +  1 

IF(J.GT.ll)  J=l 

CFDSK  J,K)=I  :iFO(  I  ) 

ICGUNT=ICGU>\T+1 

1C  =  IC«-1 

IF( ICGUNT.EQ 

CONTINUE 

GC  TO  2 

IFLAG1=1 

IFLAG=1 

IF(( IFLAG.EQ 

WRITE (6,104) 

FORMAT*  ■  ■  ,<  JOB 
116, IX, 'NO.PTS=« 

FLAG=1 

RETURN 

WRITE(6,105)  IC 


TO  15 


14 

16 

50 

51 

104 


121 )  GO  TO  16 


1) . AND.i IFLAG1 .EQ 
ICOUNT, MATRIX, IC 

FAILURE' ,1X, •COUNT=», 

16) 


1)  )  GO  TO  5  2 


16,  IX, »MATRIX  = 


=  i 


FCPMAT (« 

RETURN 
END 
AN 


','JOB  COMPLETE' , IX, • POINTS  ST0RED=',I6) 


SUBROUTINE  INTERP 
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SUBROUTINE  INTERP(ALFAC,FLAP,COEF,M) 
C C MMON  CO EFMX , CO  EFCX , C  OEFOX 

DIMENSION  CCEFMX(3,33) ,COEFCX( 3 , 33 ) , COE FOX < 3, 33 ) 
DATA  CONDEG/57. 295773/ 
ALFA=ALFAC*CONDEG 

IF(ALFA.GT.IO.O)  WRITE(6,100)ALFAC,M 
IF(ALFA.GT.IO.O)  ALFA =10.0 

FORMAT (  »0» ,' ERROR.  ALFA    EXCEEDED    L IMITS . ALF A=» ,E14. 8, 
12X,'M=' ,11) 

WRITE  (6,  100)     ALFACM 
ALFA=-6.0 


IF(ALFA.LT.-6.0) 

IF(ALFA.LT.-6.0) 

ALFHA=2.0*ALFA 

IFULPHA.LT.O.O  ) 

IF( ALPHA. GT. 0.0) 

IA=13 

WTA=0.0 

GC    TO    3 

IA1=-ALPHA 

IA=13-IA1 

£LFHA1=-IA1 

WTA=ALPHA~ALPriAl 

GO  TO  3 

IA1=ALPHA 

I A=l 3+ I Al 

ALFHA1=IA1 

WTA=ALPHA-ALPHA1 

IF(FLAP.LT.-5.0) 

IF(FLAP.LT.-5.0) 

IF(FLAP.GT.5.0) 

IFCFLAP.GT.5.0) 


GO 
GO 


TO 
TO 


WRITE(6,101)  FLAP 

FLAP=-5.0 
WRITE(6,101)  FLAP 
FLAP=5.0 
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101    FORMATC  ,0« ,»ERROR.FLAP    EXCEEDED    LI  MI TS. FLAP=« ,E14 . 8) 
IF(FLAP)4,5,6 

4  IB=1 

WTB=-< ( 5.0+FLAPJ/5.0) 
GO    TO    7 

5  IB  =  2 
WTB=0.0 
GO    TO    7 

6  IB  =  3 
WTB=(5.0-FLAP  )/5.0 

7  I F(WTA) 8,9,10 

8  IA2=IA-1 
GC    TO    11 

9  IA2=IA 
GC    TO    11 

10  IA2=IA+1 

11  IF(WTB) 12, 13,12 

12  IB2=2 

GO    TO    14 

13  162=13 

14  GC    T0<15,16,17)  ,  M 

15  C0EF1=C0EFMXU3,  IA) 
C0EF2=C0EFMX( IB2, IA) 
GC    TO    18 

16  CCEF1=CCEFCX(IB, IA) 
CCEF2  =  C0EFCXU B2 ,  IA) 
GC    TO    18 

17  C0EF1=C0EF0X( IB, IA) 
CCEF2=CCEF0X( IB2, IA) 

18  C0EF3=(C0EF2-C0EF1 ) * ABS( WTB ) +COEF 1 
GO    TO    ( 19,20,21 )  ,M 

19  CCEF1P=C0EFMX(  IB,  IA2) 
C0EF2P=CGcFMX< IB2,  IA2) 
GC    TO    2  2 

20  C0EF1P=C0EFCX< 13, I A2) 
CCEF2P  =  C0EFCX( 162,  IA2) 
GC    TO    22 

21  CGEF1P  =  C0EF0XU3,IA2) 
C0EF-2P--C0EF0X(  18  2,  1A2) 

2  2       CCEF3P=<CCEF2P-C0EF1P»*ABS(  WT3H-CCEF1P 

CCEF=(C0EF3P-C0EF3)*ABS(WTA)+C0EF3 

RETURN 

END 
FORTRAN 

C  SUBROUTINE     INTRP1 

C 

c 
c 

SUBRGUT  INE  INTRP1 (BETA1 , SOB ,COEf S, H, THET A, THETA1  1 
CCMMON  JUNK, CFDS 1 , CADS  1 

DIMENSION  JUNK (297) , CFDS1 ( 1 1 , 1 1 )  , CAD SI ( 1 1 , 11 ) 
DATA  CONDEG/57. 295773/ 
SUB=SOB 

BETA=BETA1*CGNDEG 

IF(BETA.GT.15.0)  WRITE ( to, IOC)  BETA1,M 
IF(EETA.GT.15.0)  BETA=15.0 
100  FORMAT! ■O1 ,« ERROR. BETA  EXCEEDS  L I MITS. BET A= • , E14. 8, 2X, 
l'f'=«  ,  II  ) 
IF(BETA.LT.-15.0)  WRITE(6,100>  BETA1,M 
IF(BETA.LT.-15.0 )  BETA=-15.0 
IF(BETA.LT.O.O)  GC  TO  1 
IF(BETA.GT.O.O)  GO  TO  2 
IA=6 
WTA=0.0 
GC  TO  3 

1  IA1=-BETA 
IA=6-IA1 
BETAG=-IA1 
HTA=BETA-BETAG 
GC  TO  3 

2  IA1=BETA 
IA=6+IA1 
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//L 
//L 
//G 
//G 

0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 


BE 
WT 
IF 
IF 
IF 
IS 
WT 
GO 
IS 
IS 
St 
WT 
IF 
IA 
GC 
IA 
GC 
IA 
IF 
IS 
GC 
IS 
IF 
CO 

cc 

GO 
CC 
CC 
CO 
IF 
CC 
CC 
GC 
CC 
CC 
CC 
CO 
TH 
PE 
EN 

.SYSU 

.SYSL 

.FT06 

.AUXD 

18 

24 

19 

17 

24 

135 

17 

23 

16 

16 

22 

165 

16 

22 

17 

15 

21 

17 

15 

20 

17 

14 

20 

175 

14 

19 

18 

13 


10 
11 


12 

15 


16 
18 
19 


20 
22 


TAG=IA1 

A=BETA- 

(SU8.LT 

(SUb.GT 

(SUB.GT 

=  1 

S  =  0.0 

TO  7 
1  =  SUB 
=  I S 1  + 1 
B1=IS1 
S=SU8-S 
(WTA) 8, 
2=IA-1 

TO  11 
2=IA 

TO  11 
2=IA+1 
(WTS.GT 
2=1  S 

TO  15 
2=IS+1 
( M  .  E  C  .  2 
EF1=CAD 
EF2=CAD 

TO  18 
EF1=CFD 
EF2=CFD 
EF3=(CC 
(M.EQ.2 
EF1P=CA 
EF2P=CA 

TO  2  2 
EF1P=CF 
EF2P=CF 
EF3P=  (C 
EFS=(CO 
ETA1=TH 
TURN 
D 

Tl  DO  S 
MOD  DO 
F001  DD 
ATA  DD 


BFTAG 

.0.0)  SUB=0.0 
.10.0)  SUB=10.0 
.0.0)  GO  TO  5 


UB1 
9,10 


,0.0)  GO  TO  1 


)  GO  TO  16 
Sl(  IStlAl 
Sit  IS2i IA) 

Sl( IS, I Al 

Sl( IS2, IA) 

EF2-C0t:Fl)*WTS  +  C0EFl 
)  GO  TO  20 
DS1 (  IS, IA2) 
DSKIS2,  IA2) 

DSK  IS,  IA2) 
DSH  IS2,  IA2) 

0EF2P-C0EF1P)*WTS+C0EF1P 
EF?P-C0EF3)*A6S<  WTAJ+C0EF3 
ETA*CQNDEG 


PACE=<  3500,  (  100, 10)) 
SPACE=(3500,(  100,10,1) i 
SPACE=(CYL, ( 3i J 
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0.016 

0.018 

0.013 

0.018 

0.019 

0.013 

0.018 

0.020 

0.013 

0.018 

0.020 

0.013 

0.019 

0.021 

0.013 

0.020 

0.022 

0.0135 

0.020 

0.023 

0.014 

0.022 

0.025 

0.014 

0.023 

0.026 

0.015 

0.024 

0.028 

0.016 

0.026 

0.030 

0.0175 

0.0275 

0.0^2 

0.020 

0.03C 

0.035 

0.021 

0.0^2 

0.0375 

0.023 

0.034 

0.040 

0.025 

0.037 

0.044 

0.028 

0.040 

0.0475 

0.030 

0.044 

0.052 

0.034 

0.047 

0.055 

0.036 

0.051 

0.059 

0.040 

0.054 

0.068 

0.044 

0.058 

0.067 

0.047 

0.063 

0.070 

0.050 

0.065 

0.074 

0.025 
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0.026 
0.027 
0.024 
0.025 
0.026 
0.023 
0.023 
0.024 
0.020 
0.021 
0.022 
0.018 
0.019 
0.021 
0.0165 
0.0185 
0.020 
0.015 
0.017 
0.013 
0.014 
0.015 
0.0175 
0.013 
0.0145 
0.013 
0.013 
0.014 
0.017 
0.011 
0.013 
0.018 
0.010 
0.012 
0.018 
0.008 
0.012 
0.019 
0.007 
0.012 
0.020 
0.006 
0.012 
0.021 
0.005 
0.013 
0.023 
0.005 
0.0135 
0.028 
0.005 
0.014 
0.026 
0.0055 
0.015 
0.027 
0.006 
0.016 
0.030 
0.0065 
0.0175 
0.0325 
0.007 
0.020 
0.0  35 
O.OOo 
022 
0  38 
010 
0245 
0.043 
0.012 


0 
0 

0 
0 
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0.0265 
0.047 
0.014 
0.0295 
0.054 
0.016 
0.033 
0.060 
0.018 
0.037 
0.07C 
0.021 
0.042 
0.078 
0.028 
0.048 
0.084 
0.026 
0.056 
0.093 
0.032 
0.067 
0.  10 
0.036 
0.08 
0.11 
0.021 
0.021 
0.0215 
0.02C5 
0.0205 
0.021 
0.020 
0.020 
0.02C5 
0.018 
0.018 
0.019 
0.017 
0.0175 
0.018 
0.0165 
0.017 
0.0175 
0.015 
0.0155 
0.016 
0.014 
0.0145 
0.0155 
0.0135 
0.014 
0.015 
0125 
013 
0145 
01C 
0.012 
0.015 
0.0095 
0.011 
0.015 
0.008 
0.0105 
0.016 
0.0075 
0.01C 
0.017 
0.006 
0.0105 
0.0195 
0.005 
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O.Oil 

O.C22 

0.0045 

0.012 

0.024 

0.0045 

0.013 

0.028 

0.005 

0.0145 

0.031 

0.005 

0.0165 

0.035 

0.006 

0.020 

0.041 

0.003 

0.023 

0.048 

O.O10 

0.027 

0.055 

0.012 

0.031 

0.06  4 

0.015 

0.036 

0.078 

0.019 

0.044 

0.093 

0.022 

0.051 

0.  10 

0.028 

0.060 

0.  10 

0.035 

0.075 

0.10 

0.039 

0.084 

0.10 

0.045 

0.092 

0.  10 

0.05 

0.095 

0.10 

0.056 

0.10 

0.10 

0.25 

0.30 

0.345 

0.41 

0.46 

0.52 

0.56 

0.62 

0.65 

0.70 

0.78 

0.25 

0.  295 

0.34 

0.38 

0.425 

0.47 

0.52 

0.55 
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0. 

,60 

0. 

,63 

0. 

,68 

0. 

,25 

0. 

,29 

0. 

,325 

0. 

,36 

0. 

,40 

0. 

,43 

0. 

,47 

0. 

,50 

0. 

,545 

0. 

,575 

0. 

62 

0. 

,25 

0. 

,28 

0. 

,32 

0. 

,345 

0, 

,37 

0. 

,41 

0. 

44 

0. 

,47 

0. 

,50 

0, 

,54 

0. 

,57 

0. 

,25 

0. 

,23 

0. 

,31 

0. 

,34 

0. 

,36 

0, 

,40 

0. 

,43 

0. 

.45 

0. 

,48 

0. 

,51 

0, 

,54 

0. 

,25 

0. 

,28 

0. 

,31 

0, 

,34 

0. 

,36 

0, 

,39 

0. 

,42 

0, 

,45 

0, 

,47 

0. 

,50 

0, 

,53 

0. 

,25 

0. 

,28 

0. 

,31 

0, 

,34 

0. 

,36 

0. 

,40 

0. 

43 

0. 

45 

0. 

,48 

0. 

,51 

0. 

54 

0, 

,25 

0. 

28 

0. 

,32 

0, 

345 

0. 

37 

0, 

,41 

0. 

,44 

0, 

,47 

0. 

,50 

0. 

,54 

0. 

,57 

0. 

,25 

0, 

,29 

0. 

,325 
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0.36 

0.40 

0.43 

0.47 

0.50 

0.545 

0.575 

0.62 

0.25 

0.295 

0.34 

0.38 

0.425 

0.47 

0.52 

0.55 

0.60 

0.63 

0.68 

0.25 

0.30 

0.345 

0.41 

0.46 

0.52 

0.56 

0.62 

0.65 

0.70 

0.78 

0.05 

0.11 

0.17 

0.23 

0.29 

0.35 

0.40 

0.46 

0.52 

0,58 

0.64 

0.05 

0.10 

0.15 

0.195 

0.245 

0.29 

0.34 

0.38 

0.44 

0.48 

0.54 

0.05 

0.095 

0.14 

0.  17 

0.21 

0.25 

0.28 

0.33 

0.36 

0.40 

0.45 

0.05 

0.0  8 

0.13 

0.  15 

0.18 

0.22 

0.25 

0.28 

0.32 
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0.35 

0.38 

0.05 

0.08 

0.11 

0.145 

0.  17 

0.20 

0.23 

0.26 

0.29 

0.33 

0.35 

0.05 

0.07 

0.  11 

0.14 

0.17 

0.  195 

0.22 

0.25 

0.325 

0.30 

0.33 

0.05 

0.08 

0.1  1 

0.145 

0.  17 

0.20 

0.23 

0.26 

0.29 

0.33 

0.35 

0.05 

0.08 

0.  13 

0.15 

0.18 

0.22 

0.25 

0.28 

0.32 

0.3  5 

0.38 

0.05 

0.095 

0.14 

0.17 

0.21 

0.25 

0.28 

0.33 

0.36 

0.40 

0.45 

0.05 

0.  10 

0.15 

0.  195 

0.24  5 

0.29 

0.34 

0.38 

0.44 

0.48 

0.54 

0.05 

0.  11 

0.  17 

0.23 


124 


0.29 

0.35 

0.40 

0.46 

0.52 

0.58 

0.64 

//PLCT.SYSIN  CD  * 

RESPCNSES    WITH    HELM 

RIBEIRO 

0.0  300.0 


STEP    COMMAND 


-500.0 


300.0 


5.0 


6.0 
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